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Abstract 



I The modified Kadomtsev-Petviashvili (mKP) equation is shown in this paper to be 

I decomposable into the first two sohton equations of the 2A^-coupled Chen-Lee-Liu and 

' Kaup-Newell hierarchies by respectively nonlinearizing two sets of symmetry Lax pairs. 

In these two cases, the decomposed (l-l-l)-dimensional nonlinear systems both have a 
^ ' couple of different Lax representations, which means that there are two linear systems 

' associated with the mKP equation under the same constraint between the potential 

and eigenfunctions. For each Lax representation of the decomposed (l-l-l)-dimensional 
nonlinear systems, the corresponding Darboux transformation is further constructed 
such that a series of explicit solutions of the mKP equation can be recursively generated 
with the assistance of symbolic computation. In illustration, four new families of 
solitary-wave solutions are presented and the relevant stability is analyzed. 
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1. Introduction 

According to Lax's tlieory [l'], a given NLEE is said to be integrable if it arises as the 
compatibility condition of two linear eigenvalue equations which are usually called a Lax 
pair and comprised of the spatial part and the temporal part. Although it is not an easy 
work to find the Lax pair associated with an integrable NLEE, one can relate a properly- 
chosen spectral problem to a hierarchy of soliton equations whose Lax pairs have the same 
spatial part but different temporal parts [2I . In the past several decades, some important and 
typical (l-l-l)-dimensional integrable hierarchies have been established and well understood, 
including the Ablowitz-Kaup-Newell-Segur (AKNS) jl], Wadati-Konno-Ichikawa (WKI) 0], 
Kaup-Newell (KN) jsl and Levi [gI hierarchies. Today, the Lax pair has been playing a 
considerable role in studying the integrable properties of NLEEs such as the Hamiltonian 
structures, conservation laws and symmetry classes 0]. 

For describing various complex nonlinear phenomena of our realistic world, the higher- 
dimensional NLEEs appear very attractive in many fields of physical and engineering sci- 
ences jsl. However, due to the higher space dimensions, those higher-dimensional nonlinear 
systems often exhibit more intricate properties (e.g., the integrability of NLEEs in 3+1 
dimensions is still not a well-solved problem ^) and admit more abundant soliton struc- 
tures It is mentioned that by dimensional splitting the higher- dimensional problem can 
be reduced to several lower- dimensional ones which are easier to treat with the available 
tools [ill, [12]. In recent studies, many (2+l)-dimensional integrable NLEEs have shown to 
be relevant with some known (l+l)-dimensional soliton equations by the nonlinearization 
of their Lax pairs and adjoint Lax pairs H, H Q, E^ . Through such a decomposition, a 
submanifold of solutions for the given (2-1-1 )-dimensional integrable NLEE is obtainable by 
solving the resulting (l-l-l)-dimensional integrable systems. 

The modified Kadomtsev-Petviashvili (mKP) equation Q| , as one of the most important 
integrable NLEEs in 2+1 dimensions, has been derived in many physical applications such 
as the propagation of ion-acoustic waves in a plasma with non-isothermal electrons [l3] and 
the electromagnetic wave description in an isotropic charge-free infinite ferromagnetic thin 
film 181]. It has been found that the mKP equation is able to be decomposed into the first 
two nontrivial nonlinear systems in the Burgers hierarchy 15|], two-coupled Korteweg-de 
Vries (KdV) hierarchy 19], two-coupled Chen-Lee-Liu (CLL) hierarchy 15,Eq|, two-coupled 
KN hierarchy 2l|, |22] and some other soliton hierarchies 23|]. But to our knowledge none 
has given a systematic way of determining its all decompositions to two (l+l)-dimensional 
integrable NLEEs in the same hierarchy, which means that some new or more generalized 
integrable decompositions for the mKP equation have not been uncovered as yet. 

In Ref. 15|, the authors have pointed out that the following mKP equation 
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can be constrained into the two-coupled CLL and high-order CLL systems by imposing 
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the nonlinearization on both the associated Lax pair and another auxiliary Lax pair. The 
present paper is intended to make a further investigation on Eqn. f lLip by proposing two 
types of generalized integrable decompositions which respectively reduce Eqn. f lLip to the 
first two nontrivial members in the 2A^-coupled CLL and KN hierarchies. On this basis, 
our next concern is to derive the Lax representations of the decomposed (l+l)-dimensional 
nonlinear systems and construct their respective Darboux transformations by which some 
new solitary- wave solutions are expected to be revealed for Eqn. (11.11) . 



2. Proposal of generalized integrable decompositions 

As indicated in Refs. E^ . Eqn. (11. ip is associated with the following linear system 

Uy = Liu, Li=dl — 2qdx, (2.1a) 

Ut = M,u, Mi=dl-3qdl + ^{q^-q,-d-\y) d^, (2.1b) 

from which we follow the definition of the adjoint of a differential operator (for a differential 
operator in the form Q = T^akd^, its adjoint form is taken as Q* = T,{—dx)'^ at, where the 
asterisk denotes the adjoint operator [l^) and obtain other three linear systems respectively 
with respect to v = (ux)*, m = u* and p = v*, as follows: 

Vy = L2V, L2 = -dl-2qdx, (2.2a) 

Vt = M2V, M2 = dl + 3qdl + ^ {q^ + - d'^qy) d^, (2.2b) 

ruy = — L* m, L\ = d1 + 2qdx + 2 g^., (2.3a) 

mi = -Ml m, M* = - 93 - 3 g 92 - ^ (g2 + 3 g, - d-\y) 

3 

+ l^ily -'^QQx - Qxx) , (2.3b) 

Py = -Llp, L; = -dl + 2qdx + 2qx, (2.4a) 

pt = -M;p, M; = -dl + ?>qdl-^ (g2 - 3g,. - d-^qy) 

3 

+ -(gj,-2gg,. + g,,), (2.4b) 

where System (12.20 has also been given in Ref. 15(], while Systems (12. 3p and (12. 4p are 
exhibited here for the first time. By direct calculations, we find that the compatibility 
conditions Vyt = Vty, rriyt = rrity and Pyt = Pty all give rise to Eqn. (11.10 . which suggests that 
Systems (I2.2l) - (l2.4p are other three different Lax pairs of Eqn. (II. ip . 

Assuming Uj and Vj (j = 1,2, ...,A^) respectively satisfy Systems (12. ip and (12.20 . we 
introduce the following potential constraint 13 1 
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into Systems (12.11) and (12. 2p . and obtain the following 2A^-coupled CLL system [2J] 
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j,xx UkVkUj^x = 0, Vj^y + Vj^xx ~y^^VkUk Vj,x = 0, (j = 1,2, . . . , A^), (2.6) 



k=l 



k=l 



and its high-order generalization 
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^UkVk\ + 2 ^ Uk,x Vk 
vfc=l / k=l 



'^j,x 0, 



(j = l,2,...,iV), (2.7a) 



3 ^ 3 

'^j,t Vj,xxx ~l~ 77 / J '^k Uk VjjXX T 



k=l 
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N 



Vk Uk 

,fc=i 



2 ^ ^tfc 
fc=i 



= 0, 

(j = l,2,...,iV). 



(2.7b) 



Similarly, if we constrain the potential as 



1 ^ 



{21 



where ruj and (j = 1, 2, . . . , N) satisfy Systems (12.31) and (12.4 1 ). respectively, then Eqns. (I2.3al) 



and (I2.4ap are nonlinearized into the 2A^-coupled KN system |5|, 
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'^j,y + f^j,xx - mkPkrrij ] = 0, pj^y - p^^.^. - PkrrikPj \ = 0, 



.A;=l 



.fc=l 



(j = l,2,...,iV), (2.9) 



and Eqns. (12.3bp and (12.4bp become 



ruj^t - rrij, 



N 



N 



N 



^mkPk] rrij -^mtPk rrij^^ - ^ nik^x Pk mj 



,k=l 



k=l 



k=l 



(j = l,2,...,iV), (2.10a) 



Pj,t - Pj,. 



N 



N 



N 



^Pkmk\ Pj + ^ Pk ruk Pj^x + Pk,x rrik Pj 



k=l 



k=l 



(j = l,2,...,iV), (2.10b) 



which is a generalized high-order version of System (12. 9p . Without any difficulty, one can 
check that the above two constrained potentials qj and qjj both satisfy Eqn. (II. ip exactly. 
Thus, we have got two generalized integrable decompositions for Eqn. (II. ID . as follows: 
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Decomposition I: If {uj, vj) (j = 1, 2, . . . , A^) is a compatible solution of Systems (12.60 
and (12.71) . then the function qj determined by Expression (12. 5p solves the mKP equation (II. ip . 

Decomposition II: If {mj,pj) (j = 1, 2, . . . , A^) is a compatible solution of Systems (12.90 
and (12.1 op . then the function qjj determined by Expression (12. 8p solves the mKP equa- 
tion ffLT]) . 



Note that the two proposed decompositions in Refs. [l5|, l20| and Refs. [2l|, |22|, respec- 



tively, correspond to the special cases of Decompositions I and II when = 1. Moreover, 



the authors in Refs. 2l|, |22| have also not shown the association of the decomposition there 



with the Lax pairs for the mKP equation. 

Under the potential constraints (12. 5p and ( 12.80 . we can gain much information about the 
mKP equation (11.10 from Systems (I2.6p - (I2.70 and (I2.9p - (l2.10p by means of various known 
effective approaches. In the following, to explore more unrevealed solutions (especially the 
solitary-wave solutions) of Eqn. (11.10 . we will deal with the decomposed (H-l)-dimensional 
nonlinear systems by employing the Darboux transformation method which has been proved 
to be an excellent technique for analytically studying integrable NLEEs and soliton prob- 
lems [23] in that it gives the general procedure to recursively generate a series of explicit 
solutions including the multi-soliton solutions from an initial solution 26|]. Once the Darboux 
transformation for a given NLEE is constructed, one only needs to solve a linear differential 
system (i.e., the Lax pair with an initial potential) and perform tedious but not complicated 
algebraic operations 27j. Hereby, an obvious advantage of the Darboux transformation 
lies in that the iterative algorithm is purely algebraic and very computerizable by virtue of 
symbolic computation 



28 
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32]. 



3. Lax representations and Darboux transformations of Systems (12.61) and (12.71) 

It is possible that an integrable nonlinear system could be associated with several linear 
spectral problems, which might lead to different Darboux transformations. We consider the 
(A^-|-l) X (A^+1) linear eigenvalue problem and interestingly find that Systems (12.60 and (12.70 
admit two different kinds of Lax representations, in which the first one is of the form 



= 1/(1) ^ 



(1) 



2ti/(i) + ;^w(i) 



where A is the eigenvalue parameter, ^ 



(^"1,^2 



(3.1a) 
(3.1b) 
(3.1c) 



the vector transpose) is the vector eigenf unction, the matrices ul'^\ V^^^ and 
A; = 0, 1, 2; / = 0, 1, 2, 3) are expressible in the form 



(1) 



(1) 



(1) 




(the superscript T denotes 

= 0,1; 

(3.2) 
(3.3) 
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\VA+\V {U^V - UV,) - V^UV + v,, -\ vuvu + v^u - vu. 



1 



3 



2 

-^AU+iUV^-U^V)U-2UVU^-2U^^ 
-i VAU + V.UVU - VUVU, + i V{UV, - U,V) U + V,U, - V,,U - VU, 



,(3.6) 



with / as the NxN identity matrix, A = UVUV, U = {ui, U2, ■ ■ ■ , un), V = (fi, f2, • • • , vn) 
and Oi = (0, 0, . . . , 0)^. Based on the matrix-form inverse scattering formulation in Ref. 29], 
the second Lax representation of Systems fl2.6p and (12.71) is presented as follows: 

= ?7(2) ilr = [x^ f/(2) + X f/f ) + f/f (3.7a) 

q^, = = [x^ wi'^ + X' lyf ) + lyf ^ + + \' W^^ + A + W^] ^, (3.7c) 

with 

lyf ) = 4 \ ^ = -2 ^2^^^ W^f ) = -2 1^3^^^ (3.9) 





' \0i VUj' ^ yvt/V-K ^1' V ; 

V;^'^=f° 1 (3-12) 

<)=H^ + ^^^- ^^-'^ I, (3.13) 
= 



(3.14) 



\AU + l {U^V - f/K.) + UVU, + 

j - i F([/\4 - f/,\/) - + O2 

lyf ) = i X 
^ 2 

Of 

-i VAU - VU,, - V,,U + V,U, - VUVU, + V,UVU + i V{UV, - U,V) u^ 

where O2 is the NxN zero matrix. A, I, U, V and Oi have been defined as above. Here, 
it is easy to verify that Systems (12. 6p and (12.71) can be respectively derived from the zero- 



,(3.15) 
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curvature conditions f/^*^ - vi'^ + [U'^^ , V^^] = and U^'^ - W^'^ + \U^\ = (z = 1, 2), 

where the brackets represent the usual matrix commutator. 

It is known that the Darboux transformation is actually a gauge transformation which 



relates two different solutions of the same linear system [25||. Starting from System (13.11) . we 
construct the first Darboux transformation for Systems (12. 6p and (12. 7p in the form 



4(1) Aixr 



0(1) 0(1) / ' 

D3 D4 



(3.16) 



where 



(1) _ ^(1) ^(1) 
"11 ' ^2 



(1) 



.(1) q(1) 
'11 ) ^2 



.(1) 
'12 ) 



(1) 

12 ' • • 

Ji) ^ c 



x(i) ^ ^(1) _ 



5(1) 

■^21 ' 



^(1) 



(1) 
4 



(1) ^ (,(1) 



21 ' 



.(1) 
'Af+1 



(1) 
4 



2<i,fc<Af+l' 



and sJP (1 < i,A; < A^+1) are all the functions of x, ?/ and t to be determined, If' is required 
to satisfy System dSl]) with f/^^), l/^^^ and W^^^ replaced respectively by f/^^^, V^^^ and W^^^) 
in which the old potentials (uj, f j) are transformed into new ones (wj, Vj) {j = 1,2,..., A^). 



2<i,fc<Af+l' 



ik 



From the knowledge of the Darboux transformation, we know that the matrices A^^\ S^^\ 

0, 1, 2; / = 0, 1, 2, 3) must satisfy the following equations: 



ul'\ v^'^ ^nd w; 



(1) 



0,1; k 



Zi(i)f/«_?7^^)Z\« = 0, 

zid)^/!) _ zi(i)5(i) vi''> + Vo(')zi(i)s(i) - v/i)z\(i) = 0, 

zi(i) + zi(i)v« _ ^(1)^(1) yw + f a)zi(i)5(i) _ yw^(i) ^ 0, 

^(1)^(1) + /id) 5(1) + ^(1)5(1) yd) _ yW^WsW = 0, 



A 
A 



+ 



(1) 



zi(i)5(i) pyd) + ^-(1)^(1)^(1) _ ^(1)^(1) ^ 0, 



Zi(i)wd) _ ^(1)^(1) ^(1) + wd)^(i)^(i) _ ^(1)^(1) ^ 



5(1) + zi(i) cd) + zi(i)5(i) - #(i)z\(i)5(i) = 0. 



(3.17a) 
(3.17b) 
(3.17c) 
(3.17d) 
(3.17e) 
(3.17f) 
(3.17g) 
(3.17h) 
(3.17i) 



By Eqns. (I3.17al) and (]3.17bp . we can directly compute out: 



^(1) 
^2 

^(1) 



(1) 
4,x 



oL A 



= 02 

UA 



(1) 



A 



(1) 

3,a; 



(1) 



(3.18) 
(3.19) 

(3.20) 



with 



U={5'^U-25't^S'^'){A'l')-\ V={2A'^^ + Afv)/5^l (3.21) 
Then, using the above results, it can be found that Eqn. fl3.17dp is lead to be satisfied auto- 



matically, while Eqns. fl3.17el) . p.l7gP and (I3.17hp yield the following constraint conditions 
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on 5S>, A^^^ and A^^ as 

Zig = O2, A['} = O2, (3.22) 

+ 2 5^'^USP -A6['^Si'^S^^\ (3.23) 
+ SS>UV.S? - S['^U,VS['^ - \ 5['^AS?^ + UVUA^i^si'^ + 2 U,A^^^S['^ 

+ ^ IJVjjA^^'^ - ilVU^Af^ - ^ UjrilAf^ - ^ AIJA^^'^ + 2 f/,ZiJ^S'^'\ (3.24) 

+ 2 zi^^^c/^w + zi?V[/4^) + 4 zii^^^f ^^i^) + 2 zi«5(^V5{^) 
+ 2 zi«5« v^;^) - 1 vu^A^^ + i v;?7zi(^) - \ VUVUA^^ 

+ ziJV.^w - \ A^pvuvs? - /i^c/y^w - 2 zi?V(5f 

+ ZiJ^^S^f/y - 2 Zi«ysf + zlWs(^)[/\/ - ^^S^^V^, (3.25) 
Zig = Zl«5i%.l^ - A^^S^UV^ + \ A^^S?A - A^^S^^^UV. + A^^si'^U.V 

+ \ vfjvfjAfsf - Afs^^V^UV - \ Afs^^^VUV^ + \ AfsfvU^f 

+ zi«5?V.. + I A^^S'^^VA + Zii^^^f + Vtj.A^^^si'^ - V;^7zi«5« 

+ i - Af^si'^UV - I A^^S^l^VUV. + i A^^S^Wu^ 

+ ^ vu%uA^^^ - ^ vtJvu.A^^'^ + ^ vMuA? - \ vijMa? 

+ VU^A'i^Si^^ - I vAuA'i^ - V^IJ A^l^ S^i'^ + i A^pS^^^A, (3.26) 
8 2 

and three redundant equations: 

5^iuVS^^^ - \ 5^^UVU - 2 5ff5«V^4'^ + S^^S^WU + 4 5ff ^^^^^ 

+ i UVUA^l^ - 2 41^^751'^ + il^A^P + 2 ^il^^J'^t/ - = 0, (3.27) 
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+ 2 zil'^sf V, + 2 z\^'^4'^K - 2 Vf/Zi^'^sf ^ + 5ff vf/y^i'^ - 2 (5{l^v;^f ^ 
+ zii'V,, + vu^A'i^ - v^uA'i^ + ^ vuvuA'i^ + v;f/vz\f ^ _ <5(i)v;^ 

+ \ 5iiVlJ% - ^ 6^^^VUj/ - A'^I^sI^VUV - 2 VJJA'^^^s'i^ = 0, (3.28) 
i 5SUV^U - 5S>UVU^ - i 5[l^U^VU - UVUA'i^S^^^ - 2 U.,A^^S^^^ - 5ff f/V^S^'^ 

+ UVU^aI^ + i f/^V^f/Zii^^ + sl^^Si^'^UVU + i Sj^Vf/\/f/ + 2 5ff ^t/^ 
+ 5ff S-^'V^y, + i 5ff i^'^ - '^ff ^x.: + U^^^4^ - \ S^AU + ^ AUA^l^ = 0, (3.29) 
where A = UVUV. 

The removal of u[^\ ¥2^^ and W^^^ respectively from Eqns. fl3.17cl) . fl3.17fj) and fl3.17ip 
by virtue of other six equations in System (13.171) gives 

5« = [f/«5«,5«] + (3.30a) 
= [yo(^)(5«)',5«] + + [Vi'\s'-'^, (3.30b) 

^« = [l^«(5«)^S(i)] + [W^f)(5«)',5W] + + [W^'\S^'^]. (3.30c) 

Here, after substitution of Eqns. (I3.19p - (l3.2ip . Eqns. (I3.27p - fl3.29p are proved to be comple- 
tely covered by Eqn. (I3.30ap . In other words, Eqns. (I3.27p - (l3.29p a re all satisfied identically 
with the identity of Eqn. (I3.30ap . In a similar way like in Ref. [25|, we take S^^^ as 

^(1) =^(i)A(i)(i7(i))-\ (3.31) 

with 

if« = (/.«, hS\ . . . , h^i\,), A« = diag(A«, A«, . . . , A« 0, (3.32) 

where h^f!'^ = (/^i^ , /i2fc^ • • • ' ^iv+i fc)^ corresponds to the column solution of System ( 13. ip 
with A = A^^^ (Af ^ ^ X'j}^ when i k; 1 < i, k < N + l) , namely, 

i/CD = [/(^)ij(i)A(i) + f/f (3.33a) 
i/W = V«i7W(A«)' + y«i^WA« + (3.33b) 
i7« = l^«i7(i)(A«)' + l^(^)i7«(A«)' + l^i^)if«A« + lyi^^i/W. (3.33c) 

By straightforward substitution of Expression fl3.3ip with Eqns. f l3.32p and f l3.33p . one can 
easily verify the identity of Eqns. (I3.17cp . (I3.17fp and (]3.17ip (details ignored). Thus, col- 
lecting what have been obtained above, we get the following: 
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Darboux transformation A: The matrices U^^\ V^^^ and W^^^ liave tlie same forms as 
[/(i)^ ^(1) and W^^^ under tlie linear transformation (13.161) with Eqns. f l3.31l) -( l3.33l) . where 
A^2^ = of, A^l^ is an arbitrary invertible constant matrix, and A^^^ are determined by 
Eqns. f l3.19p . (13.201) . (I3.23p - p.26p . and the relationship between the old and new potentials 
is constructed by Transformation (13.211) . 

The second Darboux transformation of Systems (12. 6p and (12. 7p based on System (13 .Tp is 
assumed as the following form (3o| 



with 




(3.35) 




where ^ ^ll ; "2 ^ V^'l^ ' ■ ■ ■ ' ? "a ^ V"21 ? • • • 5 "W+l,!/ ^ \^ik J 2<i k<N+l^ 

M2) _ (2) c^{2) _ (2) ^ ^(2) _ / (2) (2) ^T (2) _ / (2)n ' .(2) 

and sff} {I <i,k < N + 1) are all the functions of x, y and t to be determined by a set of 
equations, as below: 

Zi{2)f/(2)_f;(2)^(2)^Q^ ^3^3g^) 

^(2) _ ^(2)^(2)^(2) ^ ^(2)^(2)^(2) _ ^(2)^(2) ^ ^g^gg^^ 
^(2) + ^(2) f;(2) _ ^(2)^(2) ^^(2) _ ^(2)^(2) ^^(2) ^ f;(2) ^(2) ^(2) 

+ f/f)z\(2)5f Zi(2) = 0, (3.36c) 

^(2)^(2) ^ ^(2)^(2)^^ ^ ^(2)^(2) ^(2) ^ ^(2)^(2) ^^(2) 

- t/f )zi(^)4^) - f/f Z\(2)42) = 0, (3.36d) 
Z\(2)42) + z\(2)42)^. + zl(2)4) t/f - t/f Z\(2)4) = 0, (3.36e) 

^(2)^.^(2) _ ^(2)^(2) ^^(2) _ ^(2)^(2) ^^(2) _^ p.j2) ^ (2) ^(2) 

+ - = 0, (3.36f) 

^(2)^.^(2) _ ^(2)^(2) y(2) _ ^(2)^(2) ^^(2) _^ ^^(2)^(2)^(2) 

+ \/i^)zi(2)5f -yf Zi(^) = 0, (3.36g) 

+ Zi(2)5f -Ff Zi(2) = 0, (3.36h) 
Af^sf + Z\(2)4) + A^'^sf + Zi(2)5f 

- Zi(2)5f '- \/f )Z\(2)4) = 0, (3.36i) 
AfS^^^ + Zi^^)^) + Zi(2)5f - Zi^^)^) = 0, (3.36j) 



y \ ' I ,y 

Zi(2)^^(2) _ ^(2)^(2) ^(2) _ ^(2)^(2) ^(2) ^ ^(2)^(2)^(2) 
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+ Tyf Z\(2)4') - H/f )^(2) ^ ^3^3gl^) 

Zi(^)W^f - Zi(2)42) ^(2) _ ^(2)^(2) ^(2) ^ ^(2)^(2)^(2) 

+ lyf Zi(2)5f - lyf Zi(2) = 0, (3.361) 

+ lyf Z\(2)4') - H^f Z\(2) = 0, (3.36m) 

^(2)^(2) ^ ^(2)^(2) ^ ^(2)^(2) ^(2) ^ ^(2)^(2) ^(2) 

- lyf Zi(2)4^) - lyf Zi(2)5f = 0, (3.36n) 

^(2)^(2) ^ ^(2)^(2) ^ ^(2)^(2) ^(2) _ ^(2)^(2)^(2) ^ ^g^gg^) 

in which lj'f'\ vj:^'' and iv/^'' are required to have the same forms as uf'\ V^j^^ and VT/^^ 
(0 < z < 2; < < 4; < / < 6) except that {uj,Vj) are replaced respectively with {uj,Vj) 
(j = l,2,...,iV). 

From Eqns. fl3.36al) and (]3.36bp . we have the following results: 

aP = 01 Af = 0^, (3.37) 

?7= (5i?f/-2 5i?4^))(Z\f)-\ V={Afv + 2Afs^^^)/5fl (3.38) 

which are further substituted into Eqns. fl3.36cl) . fl3.36fj) -f r3.36h|) and (3.36k)-(3.36m), deter- 

f2^ (2\ 

mining that 5\-{ and A\ should satisfy the conditions as 

Zi2 = 02, Zi2 = 02, 42=02, (3.39) 

^fi = S?^ V - 5i? USf^ + 2 <5f,) ^ , (3.40) 

UVlIAfsf^ - U^Afsf\ (3.41) 

+ \ UV^UAfsf^ -\a tjAfsf^ - IJVU^Afsf^ + \ 5f^ASf~^ 

+ df^iJ^VSf^ - \ IJ^Afsf - b^^llJV^Sf - lJ..Afst\ (3.42) 



with A = UVUV, and yielding the following redundant equations: 

^ 5^^^UVU + (5i?f/, - ^ UVUAf^ - IJ^Af - UVdf^sf^ + 2 df^sf^U 

-2 UAfsf^ - df^Sf^VU = 0, (3.43) 

-2 VSf ^ + Vt/Z\f = 0, (3.44) 



11 



- ^Ft/ZiJ^' - 2 ?7^Z\f ^ = 0, (3.45) 

+ 5i?v;f/\/ + ^ df^vijv^ - i 5i?v^[/,r + i z\f 5f )a - z\f 5f f/K 

+ Vt/x-Zif ^ - V^ilAf + ^ F^V^t/Zi^^^ - VUVAf + i yf/V^?7Z\f 

+ VU^Afsf - VjjAfsf - b^-:^VlJVSf + 2 V;Z\f ^ = 0, (3.46) 



where the above four equations, by substitution of Eqns. (I3.38l) - fl3.41l) . can be completely 
covered by Eqns. (13.471) and (13.481) as below. 



From Eqns. (ISdl) - (l336eD . (ICTiD - dCTjD and (3.36n)-(3.36o), we remove f/f\ Jjf^ 



'2 ) 



^(2)^ j/(2)^ W^5^-* and Wg^"* by use of other equations in System (13.361) . and obtain 

+ [t/f4^4^)], (3.47) 

+ [V^)^'^ (4V + ^/'^ (4V + + , (3.49) 

4i = [^^'^ (4V + v^/'^ (4y + v^2^'^4^ + ^i'^ 4^] s? + [T/i^), 4)] 
+[Vo^^\4)]4)(4))V [K,(^\4)]4)4)4) + [vf\4)]4) 
+[v^'\s?] (4))'4^ + [vj^\4)] (4))'+ [vf), 4^4)4) 
+[\//^^4)](4))^ (3.50) 

4] = [w^f ^ + w^i'^4^ 4^] + [w^f ^ 4^] + [wi'\ 4^] 4^ + ^i4^ 

+ Xi(4y + X24^4\ (3.51) 

41= [lyf ,4T + [w^f ,414^ + [w^f ,414^+^i(4^)' + ^2(4^)^ (3.52) 
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where Xi and X2 are defined as 

X, = [1^(^)4^)4^), 4^)] + [i^(^)(4))'4),4)] + [<)4)(4))^4)] 

+[1^(^)4)4^4)] + [<n4'0''4'^] + [<^4^4^4\4^] 
+[w^(^)(4y,4)] + (4y , 4)] + [H^f 4)4^,4)] 
+[1^^4)4), 4)] + [lyf (4))^4)] + [w^f (4))^4)] 
4^,4)] + [lyf ,4)] + [lyf 4^,4)] + [lyf ,4)] 
+[iyf)4),4)] + [lyf (4))',4)] + [lyf 4^,4)], (3.53) 
X2 = [<)4)4),4)] + [<)(4))^4)] + [<)4)4),4)] 
+[iy(^)(4V,4)] + [<)4),4)] + K)(4))^4^] 
+[w^f)4),4)] + [H^f 4\4)] + [i^f 4^,4)] 

+ [l^f,4)] + [l^f,4)]. (3.54) 

In order to make Eqns. fl3.47l) - fl3.52l) satisfied, it is sufficient to require the following relation 
(proof omitted for brevity): 

4)^^(2)^(2) (^(2)) -1 ^ ^(2) ^ i7(2)(A(2))2(/7(2))-\ (3.55) 

with 

H^^^ = {hf\ h?\ . . . , h^l,), A(^) = diag(Af\ Af, . . . , X'l,), (3.56) 

where /i^f) = (/ii^), /i2A;^ • • • ' ^tv+i fc)"^ column solution of System fl3.7p with A = A^^^ 

(Af ) ^ Af ) when i ^ fc; 1 < i, A; < + 1), namely, 

i7(2) = f/(2)i7(2) (A(2))2 + t/f )ff(2)A(2) + [/f if (2), (3.57a) 

+ V3(')i/(2)A(2) + ^^(2)^(2)^ (3_57b) 
i/f) = <)/7(2)(A(2))6 + |y(2)/7(2)(A(2))^ + W^f i7(2)(A(2))^ + lyf 1/(2) (a(2))^ 

+ lyf i7(2)(A(^))' + lyf ij(2)A(2) + PFf (3.57c) 

- (2) ~ (2) ~ (2) 

Darboux transformation B: The matrices , and Wi have the same forms as 
Uj;'^\ 14^^) and IV/^) (0 < i < 2; < A; < 4; < / < 6) under the linear transformation fIXMD 
with Eqns. (l335|) - fl337D . where Z\5^) = Of, Z\f ^ = d, Z\f ^ is an arbitrary invertible constant 
matrix, is determined by Eqns. fl3.40l) - fl3.42l) . and the relationship between the old and 
new potentials is constructed by Transformation fl3.38p . 

4. Lax representations and Darboux transformations of Systems (12.91) and ( 12.101) 

Likewise, we find that Systems (12.91) and (12.101) are also related to two different types of 
Lax representations, which are respectively written as 
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= f/(3) ,p=[X t/f ) + f/f )]^, (4.1a) 

= V(3) ^ = [A^ + A V/'^ + (4.1b) 

,p^ = W^^^<P= [X^ W^^^ + \^ W[^^ + A W!f^ + WP] (4. Ic) 

with 

^ _2 ^ ^^(3) ^ ^ ^^3) ^ ^^(3) ^ _2 ^^(3) ^ ^^^2) 

MP 






^(3)^/0 MPM-MA y^(3)^/0 |5M-|M,PM-|MPM, + M,^^ ^^^^ 
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^(3) _ / 1 5 + MP, - M,P -2MPM + 2M, 

^ ^ ^ PB + ^P^MP + lPMP^ + P^^ -^PMPM + PM^-P^M' ' 



and 

^, = f/(4)^= [A2f/(^) + Af/f)]^, (4.6a) 

= F (4) ^ = [A^ Vi'^ + X' V}'^ + \' V,^'^ + A V}'^] ^, (4.6b) 

= (4) CP = [X^ W + A^ Wi'^ + X' lyf + A=^ + A^ + A lyf ] ^, (4.6c) 

with 

K,^^) = 2 [/(^\ = 4 y/^) = 2 Ui'\ (4.7) 

(4) = 4 Ul'^ , W^i') = 2 V,^'^ , lyi") = 2 Vt^ , (4.8) 

- h > - f ° ""V - f"''"^ V (4-9) 

° \^0i / y ' \^-P O2J \ Oi PM)' ^ ' 

. ''^'i-''^]^ (4.10) 



PMPM + P,M - PM, 



I PM - I MPM,. - I M^PM + M,,. 



\P,MP -\PMP^-P,^ O2 ' 

T 



(4.12) 



where (p = (0i,02, • • • ,07V+i)^, M = (mi, ma, . . • ,mAr), P = (pi,p2, • • • ,P7v) , B = MP MP, 
Oi, O2 and / have the same definitions as in Section 3. The zero-curvature conditions 
U^^ -V^'h[U^'\ = and U^^ -W!^'h[U^'\ = (i = 3, 4) give rise to Systems fCT 

and fl2.10p . respectively. 
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Following the procedure in Section 3, we can also arrive at two Darboux transformations 
for Systems f l2.9p and f l2.10p . as follows: 

Darboux transformation C: The linear system (14.11) is kept invariant by the gauge 
transformation 

//i(3) /.(3)\ /c.{3) c-(3)\ 

i> = (AZi(3) - Zl(3)5(3)) ^, Zi(3) = (^^J3) ^^3) j , 5(3) = (^^;3^ ^^3^ , (4.13) 

where S^'^ = i7(3)A(3) (if(3))-\ A(3) = diag(Af \ A^'\ . . . , , H^'^ = {hf\ , /ig^i) 

with hf^ (1 <A;<A^ + 1) as the column solutions of System (14.11) for different eigenvalues 

X _ x(3) a(3) _ a{3) _ e(3) /,(3) _ /.(3) .(3) ^T (3) _ . (3)x , , 

following conditions: 



5fiJ^ = 2MZ\f - (5{i^MP^f^ + MPMAf - M^Af + ^f^^^f ^PMP 



= MAf + (4. 14) 

fsf^ - df^MPSf^ + MPMAf 

+ ^P^. + 2 MAfsf^ + 5!?5f ^MP, (4.15) 

= 5f^M^PS? - MP,5f ^ - ^ 5if PSf ^ + 2 MPMZif ^ + ^ PMZif 
+ 2 MPMAfsf^ - 2 M^Z\^^^5f ^ - - MPM^Af - - M^PMAf 
+ M^^Zi^^^ + ^ ^P + ^MP^ - 5f^sf^M^P - 2 M^Zif ^ 

+ \ Sf^Sf^PMP., + ^ P.MP + df^Sf^P^, + ^ 5;?4')PP, (4.16) 

Afl = 2 Z\f 5f - AfPSf^ + Z\f 5f P + Afsf^P, (4.17) 
Z\g = Afsf^PMP + Z\f 4^^^^ + /if^^^PMP + Afsf'P.^ + PMZif ^ 

+ zA^^^^f ^MP - df^PMPSf^ + Z\f 4^^^P + PMAfsP - bf}P^Sf\ (4.18) 

Zig = \ AfsfB + ^ AfsfB + Zif 5f )mP. - Zif 5f )m.P - ^Sf P..5f ^ 

+ Zif 5f P.. + \ AfsfPB + ^ AfsfPMP^ + ^ Zif ^^'^P.MP 

+ P,MZ\f - PM,^Afsf + ^MZif + ^PMPMZ\f 

+ '^PMPMAfsf - ^ (5i?P^MPSf ^ - Afsfu^P - PM^Afsf 

+ AfsfP,, - \ 6^{;PMP,S[''> - ^ 6't^PBS['\ (4.19) 



Afl = -2 Z\f 5f ^ - Zif P5f ^ - M, (4.20) 
Zi5!|J = -A^MPM - 2 Z\f - 2 Z\f 4^^M - Z\f ^^M - 5fi^P,Sf ^ 

+ PMAfsf - AfsfPM - sfi'PMPsi^^ + Af M,^ + PMAfsf\ (4.21) 
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with 



AfJ = - AfMPM^ - - AfBM + - AfM^PM - A^M^^ + 2 Af^s[^^M^ 
+ 2 Z\f ^M^. - 2 A^Si^^^MPM - 2 Z\f ^^^^PM + Z\f 
+Zif PM^. - Z\f ^P^.M - ^ Afsf^PMPM - Z\f 5f ^P^M 
+4M^f - ^ AfsfPMPM - PM.^Afsf^ - ^ PPSf^ 
+ ^ PMPMAfsf^ - ^ df^PMP^Sf^ - 5f^P^^sf^ - PM^c^fsf'^ 
+ P^MZif + ^ PMPMAfsf - ^ ^MPSf \ (4.22) 

M = (^i^M + 2 5i?^f ) (Zif )"\ P = [AfP + 2 Z\f (4.23) 



Darboux transformation D: The linear system 04.61) is kept invariant by the gauge 
transformation 

^ = (AMW - XA^^^sf - A^^^sf) ^, (4.24) 

with 

^(4)^(4)^(4) ^^(4) ^ ^(4) ^ ij(4)(A(4))2(i7(4))-l^ (4.26) 

Where A(^) = , . . . ,\%,) , E^') = {hf ,hf\ . . . Mti) with hf (1^< < 

+ 1) as the column solutions of System (14.61) for different eigenvalues A = X'f^, A^f* = Of, 
Alt^ = Oi, Zii^) = 5[t Af = (4'^)2<,_,<^+i obey the following conditions: 

C = 2 sSsi'^Slt^ - 5ff 5f P - 6['^MSi'\ (4.27) 
C = S?^MPS['^ - df^st'^PMP - P. - df^sf^MP 

+ M.Zif 5$'^ - MPMAfsf, (4.28) 
C = -\ ^'^sfPB - \ 5i?4V.MP - \ 5^^St'PMP. - ^it^sf P.. 

+ I MPM^Afst^ - 5^^M^PSt^ -^B MAfs^^ + ^ ^j^P^f ^ 

+ ^M.P - M,,4')4'\ (4.29) 

^% = Afst^M + AfPSf - 2 Zif (4.30) 
Zig = Z\f 5f MPM - Zif 4')m. + Afsf^PM - PMAfs^^ 

+ 5!?P.4') + 5!?PMP5f , (4.31) 
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+ Afst^P^M - Afsf^PM, + ^ 5ff PMP.S^') + ^ A^st^PMPM 

+^ <5i?P55f + ^ 5i?P.MP4') - ^ PAfPMZif 5f - P.MZif 5f 

+ PM.4'^5f +5ffp..4'\ (4.32) 

with 

M= (5i?M-25ff^f )(Zif)"\ P= (Zif)p-2 4')4'))/<^n- (4-33) 

5. New solitary- wave solutions with symbolic computation 

From the previous resuhs, we can gain a series of explicit solutions for the mKP equa- 
tion (11. ip by the following iterative procedure: 

(1) For the initial potentials (t/, V) and (M, P), solve the linear systems dSl]), dSZD, CT) 
and (14.61) with different eigenvalues A^*^ for column solutions h^j^ {1 < k < N + 1; 1 < i < A). 

(2) Work out the matrices S'^^\ Sf \ S'^\ S^''\ Sf \ and A^'^ (1 < i < 4), and yield 
the new potentials {U,V), {U,V), {M,P) and (M, P) via Transformations ([MI]), fl338|) . 

and (03]) . 

(3) Substitute (f/,^^), (t^, V") into Expression ([23]) and (M,P), (M, P) into Expres- 
sion (12. 8p . and obtain four families of solutions for the mKP equation (II. ip of the form 



gJJ) = -i f/V^ = -i (f/ - 2 5«) [2 {A^Ya^P + V] , (5.1) 

^S?^ = f/V^ = (f/ - 2 (V + 2 4^)) , (5.2) 

= -\ MP = -\ (M + 2Sf^) [2 {Afy'Af + P], (5.3) 

= -\ MP =-\{M-2 S^^) (P - 2 5^)). (5.4) 



Based on the above explained procedure, in what follows we shall combine Darboux trans- 
formations A-D with Decompositions I and 11 to construct explicit solutions of Eqn. (Il.ip by 
starting with the trivial seed solutions of Systems (I2.6p -( |2l71) and (I2.9p - (l2.10p . For instance, 
we solve the linear system (13. ip with = 1, mi = t>i = 0, A = A^.^"* {k = 1, 2; A^"'^'* 7^ X^) and 
get the following: 

= af^exp[\f^x + 2{xP)% + 4{\PYt], (5.5a) 
= aPexp[- X'-px - 2 {X^-^Yy - 4 {X'-p)% (5.5b) 
h^j = a«exp [X^x + 2 (A«)^y + 4 (A«)'t] , (5.5c) 
hS = a«exp [ - A«x - 2 (A«)^y - 4 (A«)\] , (5.5d) 



17 



where A^^^ A^^^ and (/ = 1,2,3,4) are all nonzero constants. From Eqn. fl3.31l) . the 
element 5^2'' matrix S^^'^ is expressed as 



(1) _ y-^l ~ ^2 J'hl 'h2 (r p\ 



Proceedingly, by performing symbolic manipulations on Eqns. (13.201) . (I3.25P and (13.261) with 
substitution of Expressions (I5.5ap - fl5.5dl) . the function is figured out as follows: 

with ^22'' an arbitrary nonzero constant. To this point, the first family of solitary- wave 
solutions for Eqn. (II. ip is obtained as 

~ (7f ^ + 7$'^ tanhe(^)) (7$'^ + 7?^ tanh^^i)) ' 
with 

= (A« - A<")x + 2[(A<'>)^ - (A<")> + i[(x'rf - (A^")']i, 

7l - - ^2 "3 ) 72 - "1 <^4 Ai - ^2 0^3 A2 , (O-yj 

73 — ttg -|- Ci]^ CK4 , 74 — '-^1 '-^4 ^1 ~r '-'2 "^3 ^2 ■ 



In a like manner, we can present other three families of solitary- wave solutions for Eqn. (II. ip : 

(5.10) 



27fH(Afy-(Ary]sech2e(^) 
7? ^sech 2 ^(2) + 7f ^ tanh 2 i^^) - 7^ 



with 

e(^) = -[(Af))^- (A?y]x + 2[(AfV- (ArV]y-4[(Af))^- (Ary]t, 

7^ =af^araraf [(AfV- (Af)^], 7^ = «f ^^r^af^af [(Af ^ + (A^)^] , (5.11) 

7f = [(«r«ry - (af^af )^] Af^Af, 7^ = + (af ^af Af ^A^ , 

(7^ + 7^ tanhe(3))(7f +7f tanhe(='))' ^' ^ 

with 

e^^) = -(Af) - AfV + 2[(Af))^ - (AfV], - 4[(Afy - {X?^Y]t, 

7f ) = af^af^ - afaf , 7^ = Af^ - af^f A?), (5.13) 

^(3) _ ^(3)^(3) , ^(3)^(3) (3) _ ^(3)^(3) 3^ (3) ^(3)^(3) 3^ (3) 
73 — Lt2 ttg -j- Ct]^ u;^ , 74 — ^4 ^1 ~r "2 3 2 ' 
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„ -2,r[(Ary-(A^'y]Bech2;W 
7fsech2eW+7<'"tanh2«('')-7f' 

with 

e(^) = [(aS^))^- (Ar))> + 2[(Af))^- (ArV], + 4[(A(^))^- (A?Y]t, 

= ai^)ar^ar^af[(AfV - (A^V] , 7^ = cf^ar^af [(aS^^)^ + (A?))^], (5.15) 

We recall that {i = 1,2,3,4) are four new families of solitary-wave solutions for 
Eqn. fll.ip since they are different from the conventional solitary-wave solutions in terms of 
a finite series of tanh and/or sech functions. Through the qualitative analysis, it is easy to 
see that the functions gf^ {i = 1,2,3,4) have no singularity and exhibit stable bell profiles 

for {af\a^\a^\ a?, A^, A?) G A^'^ U B^'\ where A^'^ and 8^'^ are defined by 

^« = {(af\a«,a«,a«,A«,A(^))|a«a« < 0,a«a« > 0,A«A« > 0, A« ^ A«}, 
= A«, A«)|a««« > 0,a«a« < 0, aS^U^ > 0, A^^^ ^ A«}, 

/I (2) r/ (2) (2) (2) (2) ,(2) ,(2)\| (2) (2) (2) (2) , (2) , (2) ^ ^ \ (2) / \ (2) "l 

^ = {("l ^"2 ^"3 ^"4 ^^1 ^^2 jl"! "2 "3 "4 < 0, A^ 'A^ ^ > 0, A^ ^ ^2 h 

= {{a?\af\4\af\xf\\f^)\af^a^^^4^af^ > 0,X?h^^^ < 0,X? ^ Af }, 
= {(af),af ,af ,af ,Af\Af af < Caf^af > CAf^A^ > 0, Af ^ ^ Af }, 
^(3) = {{af\af\af\ar,Xf\xr)\af'ar > Caf^af < 0, Af^Af > 0, Af ^ ^ A?)}, 
^(^) = ,Af\Ar))|a(^)a?)af < CAS^U?) > 0, A^^^ ^ A^^}, 

^(^) = {{af\4\af,at\X^^\X^^^)\a^^^a^^^afat^ > 0,Xfh^^^ < 0, Af ^ ^ A^^^}. 



6. Conclusions and discussions 

In this paper, we have shown that the mKP equation (11.11) can be reduced to the first 
two nontrivial nonlinear systems in the 2A^-coupled CLL and KN hierarchies, i.e.. Sys- 
tems ([2S])-([22]) and (12:91) - fl230D . by imposing the potential constraints ([23]) and ([23]) on 
Systems (I2.ip - (I2.2I) and (I2.3I) - (I2.4I) . respectively. Furthermore, it has been found that the 
2A^-coupled CLL and high-order CLL systems possess two different Lax representations (13.11) 
and (13.71) . and the 2A^-coupled KN and high-order KN systems also admit two different Lax 
representations (14.11) and (14.61) . For these four Lax representations, we have constructed the 
corresponding Darboux transformations by which abundant explicit solutions of Eqn. (II. ID 
can be obtained in a recursive manner. Through one-time iteration of those Darboux trans- 
formations, four new families of solitary-wave solutions have been presented and the relevant 
stability has been analyzed. Finally, we would like to discuss the following issues: 

1. As far as we know, there are usually two ways of finding the integrable decomposi- 
tions for a (2+l)-dimensional integrable NLEE: the first is to choose a proper (1+1)- 
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dimensional soliton hierarchy and relate its first two nontrivial members to the de- 
sired (2+l)-dimensional equation fiol . 20 . 23| : the second is to directly nonlinearize 
one or two Lax pairs of the (2+1) -dimensional integrable equation into two (1+1)- 
dimensional nonlinear systems [12, , [isl, [2o|. Obviously, the proposal of De- 

compositions I and II in Section 2 is based on the second decomposition method. 
In addition, it is noted that the potential constraint (12.51) or (12.81) . which originates 
from the Bargmann symmetry constraint 3J] , might also be applicable to some other 
(2+l)-dimensional integrable NLEEs with the availability of two symmetry Lax pairs 
like Systems (I2.ip - (I2.2I) and (I2.3I) - (I2.4I) . Special attention should be paid to that Sys- 
tems (12. ip and (12. 2p can be transformed to each other with [q{x,y,t),u{x,y,t)] < — *• 
[g(— X, —y, —t),v{—x, —y, —t)], so do Systems (12. 3p and (12. 4p with [q{x, y, t),m{x, y, t)] 
< — > [q{-x, -y, -t),p{-x, -y, -t)]. 

2. From the derivation of Darboux transformations A-D, we infer that for the following 
two general linear spectral problems 

^.= [XQil^ + Q^^^]^, (6.1) 
^y=[X'Q^o^ + XQ?+Q?]^, (6.2) 

where A is the eigenvalue parameter, ll^ = {ipi, 1^2, ■ ■ ■ , ^ = (</>!, (p2, ■ ■ ■ , (t^NY are 

the vector eigenfunctions, Qf^ and Q^^'* (i = 0, 1; = 0, 1, 2) are all the NxN matrices, 
the corresponding Darboux transformations can be respectively taken as 

^ = {XAi- AiS)^, (6.3) 
i> = {X'^A2- XA2S±- A2St)<P, (6.4) 

where Ai, A2, S, S± and Sj are five NxN undetermined matrices, and ^ re- 
spectively satisfy Eqns. (16. ip and (16. 2p with Q^^ and Q^^^ replaced by and Qf^ 
(i = 0, 1; /c = 0, 1, 2). By utilizing Ansatzs (16. 3p and (16. 4p . one can attempt to construct 
the Darboux transformations for many hierarchies of soliton equations. 

3. It is emphasized that the iterative algorithm of the Darboux transformation can be 
easily achieved on the computerized symbolic computation systems such as Mathemat- 
ica and Maple. Accordingly, if we choose > 1 and make the successive iteration of 
Darboux transformations A-D, many more complicated explicit solutions of Eqn. (11.10 
will be unearthed, which might be different from those previously obtained. 
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